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摘要 : 质子 重 离子 放射 冶 疗 中 ， 通 过 CT 扫描 获取 患者 每 个 体 素 点 的 相对 阻止 本 领 ， 是 治疗 
计划 系统 精确 计算 体内 剂量 分 布 的 关键 。 本 研究 通过 设置 常规 CT 不 同 管 电压 对 组 织 等 效 模 
体 进行 两 次 CT 扫描 ， 通 过 两 组 图 像 的 CT 值 计 算出 模 体 的 相对 电子 密度 和 有 效 原 子 序数 ， 
进而 计算 出 模 体 的 相对 阻止 本 领 。 研 究 结果 表明 ， 使 用 双 CT 值 预测 的 组 织 等 效 模 体 的 相对 
阻止 本 领结 果 优 于 常规 基于 单 组 图 像 获 得 的 结果 。 因 此 ， 本 文 方法 对 进一步 提高 质子 重 离子 
放疗 治疗 计划 剂量 计算 的 准确 性 具有 重要 作用 。 
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Abstract : In proton and heavy ion radiotherapy, the relative stopping power of each voxel point obtained by CT 
scan for patient is the key to accurately calculate the dose distribution in the body using treatment planning system. 
In this study, a tissue equivalent phantom was scanned twice by setting different tube voltages for a conventional 
CT scanner. The relative electron density and effective atomic number of the phantom were calculated by the CT 
values derived from the two-set images, and then the relative stopping power of the phantom was calculated. The 
results showed that the relative stopping power ofthe tissue equivalent phantom acquired with the two-set CT values 
was better than that based on the conventional single set of CT images. Thus, the method established in this study 
might play an important role in further improving the dose calculation accuracy in the treatment planning of proton 
and heavy ion radiotherapy. 
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Till 


1 前 


放射 治疗 (放疗 ) 作为 一 种 治疗 恶性 肿瘤 的 治疗 手段 ， 已 经 发 展 了 上 百年 。 质子 与 重 离 子 放疗 的 出 丙 
使 得 现代 放疗 又 迈 入 了 一 个 全 新 的 发 展 时 代目 。 质 子 重 离子 束 流 从 加 速 器 引出 后 射 入 人 体 ， 在 到 达 肿 瘤 病 


灶 前 ， 射 线 能 量 释 放 很 少 ， 但 是 到 达 病 灶 后 ， 射 线 会 瞬间 释放 大 量 能 量 ， 穿 过 肿瘤 组 织 后 射线 能 量 迅 速 陡 
肖 衰减 ， 此 过 程 形成 的 尖峰 称 为 布拉格 峰 (Bragg Peak)， 治 疗 时 可 根据 肿瘤 深度 选择 合适 的 束 流 能 量 ， 并 
宽 布 拉 格 峰 〈Spread-out Bragg Peak; SOBP) 巴 ， 使 布拉格 峰 履 盖 肿 瘤 ， 因 此 质子 重 离子 治疗 技术 既 可 以 提 


高 肿瘤 区 域 的 剂量 ， 又 可 以 降低 入 射 路 径 上 正常 组 织 的 受 照 剂 量 ， 降 低 放 射 毒 性 ， 实 现 疗 效 的 最 大 化 。 


E 


质子 重 离子 治疗 计划 系统 利用 计算 机 断层 扫描 〈CT) 作为 基础 输入 数据 ， 将 CT 值 (Hounsfield Unit, 
HU) 转换 为 粒子 的 相对 于 水 的 阻止 本 领 比 (Stopping Power Ratio» SPR) 或 者 水 等 效 长 度 (Water-equivalent 
path lengths，WEPL)， 从 而 计算 粒子 在 介质 中 的 射程 、 能 量 沉 积 及 剂量 分 布 69。 转 换 过 程 需要 借助 CT 值 
与 SPR 或 WEPL 转换 曲线 。CT 图 像 以 每 个 像素 点 转换 成 WEPL 后 便 可 以 将 在 水 中 计算 得 到 的 剂量 映射 到 
CT 图 像 上 。 以 水 中 剂量 分 布 曲 线 为 被 积 函数 以 每 个 像素 的 WEPL 为 积分 区 间 进 行 积 分 再 除 以 相应 的 像素 
WEPL 即 可 得 到 每 个 像素 获得 的 剂量 9。 转换 曲线 的 准确 性 直接 影响 人 体 组 织 、 嚣 官 的 阻止 本 领 的 准确 性 ， 
从 而 影响 剂量 分 布 %。 获取 CT 值 相对 阻止 本 领 转换 曲线 有 两 种 常用 方法 ,利用 动物 器 官 或 者 组 织 进行 直接 
测量 "， 或 是 利用 组 织 等 效 模 体 进行 测量 得 到 *™。 上 述 方法 带 来 的 射程 不 确定 性 为 3%-5%， 需 要 在 治疗 计 


划 中 考虑 误区 外 扩 来 减 小 影响 ™。 


iS 


带电 粒子 穿 过 介质 时 ， 其 能 量 阻止 本 领 可 用 Bethe-Block 公式 描述 ， 在 Schneider 等 "的 研究 中 ， 介 质 


的 SPR 可 用 (1) 式 近 似 表达 ， 舍 去 了 公式 中 的 高 阶 无 小 量 。 


2 (ore )- m 
SPR ~ po RE g 


(1) 


式 中 pe 为 介质 相对 于 水 的 电子 密度 ， B 为 带电 粒子 的 相对 论 速 比 系数 ，mec? 为 电子 的 静 能 量 ，I 为 介质 村 
料 的 平均 激发 能 量 ， 水 的 平均 激发 能 量 Iw=75eVe9。 


从 CD 式 可 以 看 出 SPR 除了 与 pe 呈正 比 关 系 外 ， 还 与 介质 的 工 值 呈 一 定 的 相关 性 ， 而 SPR 对 8 的 依 
赖 性 非常 小 ， 低 于 0.4%，， 在 本 文 研究 中 忽略 此 项 影响 。 根 据 中 国 科学 院 近 代 物 理 研究 所 研发 的 碳 离子 治 
疗 系统 使 用 能 量 范围 (120MeV/u 一 400MeV/) ,计算 时 统一 取 中 位 能 量 260MeV/u 时 的 B 值 ， 以 进一步 


减 小 8 值 对 SPR 的 影响 。 


在 Hünemohr 的 研究 中 ， 双 能 CT (DECT) 允许 提取 一 个 额外 的 参数 ， 从 而 解决 了 I 值 对 SPR 带 来 
的 不 确定 性 。 同 一 模 体 在 两 个 不 同 能 量 下 的 CT 值 允 许 定义 除了 pe 之 外 的 每 个 体 素 的 有 效 原子 序数 Zer”, 
以 更 准确 的 信息 来 区 分 具有 相似 pe 的 组 织 ， 从 而 使 SPR 计算 更 为 精确 。 


本 研究 利用 CT 电子 密度 模 体 中 9 个 已 知 参数 的 组 织 等 效 模 体 插件 ， 根 据 Hünemohr 等 在 研究 中 报 


道 的 DECT 双 参 数 提取 方法 ， 使 用 常规 单 能 CT (SECT) 设置 不 同 管 电压 对 组 织 等 效 模 体 进行 扫描 ， 提 取 
模 体 的 pe 与 Zes， 并 计算 SPR， 并 于 传统 方法 计算 的 SPR 做 对 比 ， 以 期 获 通过 此 方法 获得 更 加 精确 的 SPR 
值 。 


2 方法 与 材料 


2.1 设备 与 材料 


德国 西门 子 定位 CT 扫描 仪 ， 型 号 : SOMATOM Sensation Open, ŽE CIRS 电子 密度 模 体 ， 型 号 : 


CIRS 062M， 配 套 9 种 组 织 等 效 模 体 插件 ， 参 数 如 表 1: 


K 1: CIRS 062M 材料 参数 列表 


组 织 等 效 材料 pe Zeff 组 织 等 效 材料 pe Zeff 

吸 气 肺 0.2 7.332 肝脏 1.052 7.409 

呼 气 肺 0.496 7.382 骨 松 质 200mg/cc 1.117 10.55 

脂肪 0.949 6.37 骨 密 质 S00mg/cc 1.456 12.566 
g 

乳腺 0.976 6.766 骨 密 质 1250mg/ec 1.695 13.456 

肌肉 1.043 7.383 


2.2 模 体 扫描 
将 组 织 等 效 模 体 插件 放置 在 配套 的 卡 槽 内 ， 使 用 定位 CT 进行 扫描 ， 设 置 4 组 扫描 序列 ， 管 电压 分 别 


设置 为 80kVp、100kVp、120kVp 和 140kVp， 管 电流 统一 设置 为 320mAs， 扫描 层 厚 3mm， 图 像 重 建 层 


厚 1.5mm。 扫 描 示意 图 见 图 1 。 


1. 模 体 扫 描 示 意 医 


利用 RadiAnt DICOM Viewer 医学 影像 浏览 器 软件 的 圆 形 选取 工具 ， 测 量 4 组 影像 序列 中 各 组 织 等 效 


模 体 中 间 层 的 平均 CT 值 及 标准 偏差 ， 结 果 见 表 2。4 组 序列 被 测量 区 域 空间 位 置 完 全 一 致 。 
表 2. 组 织 等 效 模 体 平均 CT 值 及 标准 偏差 


组 织 等 效 材料 80kVp 100kVp 120kVp 140kVp 
吸 气 肺 -816.49+6.81 -817.69+6.15 -818.12+5.81 -817.19+5.68 
PAi -499.01+6.80 -500.03+3.78 -500.33+3.28 -500.53+3.29 
“D -97.53+5.99 -85.12+3.77 -76.68+3.06 -72.724+2.74 
乳腺 -47.89+5.34 -40.09+3.46 -34.33+2.86 -31.62+2.42 
LA 45.3247.16 44.59+5.18 45.16+4.26 43.57+4.01 
FE 57.38+7.56 55.83+4.55 56.09+3.57 53.94+3.43 
骨 松 质 200mg/cc 355.36+11.08 291.71+9.12 255.78+7.52 235.55+7.36 
骨 密 质 800mg/cc 1377.53413.05 1131.60-9.87 986.01+8.34 888.44+10.05 
骨 密 质 1250mg/cc 2119.86+16.26 1738.62+11.21 1510.78+10.2 1343.12+12.71 
2.3 参数 提取 方法 


在 Hünemohr 的 研究 中 ,将 人 体 材料 的 X 射线 线 衰减 系数 hs 近似 分 解 为 一 个 虚拟 的 纯 光电 吸收 材料 和 
一 个 虚拟 的 纯 康 普 顺 效应 材料 ， 


zn 
Ug ~ Apefg + Bpe z (2) 


其 中 n ABERRI 5 JT HUS IRI fee Y X 射线 能 量 近 似乎 坦 的 函数 ， 在 本 文 研究 中 ， 
可 视 为 常数 。A、B 是 与 材料 几乎 无 关 的 比例 因子 ，Z 为 材料 原子 序数 , 也 为 X 射线 能 量 。(2) 式 对 CT 的 
X 射线 连续 谱 也 适用 。 


对 于 CT 的 X 射线 连续 谱 ， 介 质 的 平均 线 误 减 系数 三 与 CT 图 像 的 HU 值 有 如 下 关系 : 
HU = (==) x 1000 (3) 


AB. ug. KP ISAS ERR. AER, KAEH FARER RAEE VG 73H « 


通过 高 CHO 低 (LO 两 档 能 量 扫描 组 织 等 效 材料 ， 有 如 下 表达 式 : 


g 


Hg, ~ Apefg, + Bpe zur 


(4) 


Hg, ~ Apefg, + Bpe zur 


结合 (3) x. 确定 参数 A、 B 后 ， 联 立 (3 )、 (4) 式 ， 可 以 求解 出 pe、 Zett» 并 将 已 知 数据 代入 ， 可 得 
出 如 下 结果 : 


pe = Cx (E 1) ao) x (£81) 
SET De 


由 (4) 式 可 知 ， 通 过 1 个 已 知 组 织 等 效 模 体 的 双 能 量 扫描 ， 就 可 以 确定 参数 C、D。 对 于 Z= 1-20 的 
范围 和 使 用 的 CT 光谱 ， 最 佳 系数 为 n=3.1"。 


在 本 文 研究 中 ， 使 用 全 部 9 种 材料 的 CT 影像 数据 ， 其 中 高 能 量 统一 采用 140kVp， 低 能 量 分 别 采用 


80kVp、100kVp， 拟 合 出 (5) 式 中 的 C、D 参数 ， 见 表 3。 


表 3. 参 数 C、D 拟 合 结果 


80/140 100/140 
C -0.8508 -1.6831 
D 5506.5149 10888.2124 


将 表 2 测量 结果 代入 5) 式 ， 计 算 pe、Zer， 并 与 表 1 的 已 知 参考 数据 做 比较 ， 结 果 见 表 4: 


表 4. 组 织 等 效 模 体 p。 Zor 提取 结果 与 偏差 


pe 计算 值 Zeff 计算 值 pe 偏差 Ze 任 偏差 

组 织 等 效 材料 80/140 100/140 80/140 100/140 80/140 100/140 80/140 100/140 

吸 气 肺 0.182+0.012 0.184+0.018 7.557+1.166 7.295+2.438 -9.00% -8.00% 3.07% -0.50% 

呼 气 肺 0.498+0.008 0.499+0.011 7.53540.376 7.506+0.497 0.40% 0.60% 2.07% 1.68% 
脂肪 0.948+0.007 0.948::0.010 6.54540.229 6.4994:0.326 -0.1196 -0.1196 2.75% 2.03% 
乳腺 0.9824-0.006 0.983+0.009 6.948+0.174 6.932+0.250 0.61% 0.72% 2.69% 2.45% 

几 肉 1.04240.010 .04240.014 7.49740.197 7.504+0.312 -0.10% -0.10% 1.54% 1.64% 

TAE 1.051+0.009 .05140.012 7.54140.196 7.54940.265 -0.10% -0.10% 1.78% 1.89% 

骨 松 质 200mg/cc 1.1340.017 .14140.025 9.66540.195 9.52740.353 1.52% 2.1596 -8.39% -9.70% 
骨 密 质 800mg/ce 1.47240.022 .47940.032 12.441+0.114 12.37540.190 1.1096 1.5896 -0.99% -1.52% 
骨 密 质 1250mgycc 1.682+0.027 .67740.039 13.57440.112 13.61040.174 -0.77% -1.06% 0.88% 1.14% 


2.4 InI 值 计算 
(1) 式 中 的 InI 值 使 用 上 述 方法 中 提取 的 Zor 进 行 估算 外 ， 其 与 Zor 的 关系 如 下 : 


= 0.125Zeff +3.379  Z,4,€85 a 


3.1 SPR iT 


根据 表 4 的 计算 结果 ， 结 合 〈1) 式 计算 出 材料 的 相对 阻止 本 领 ， 并 于 参考 值 作 比 较 ， 结 果 见 表 S. 


表 5.SPR 测量 结果 与 误差 


SPR 测量 相对 误差 
组 织 等 效 材料 SPR 参考 
80/140 100/140 80/140 100/140 
吸 气 肺 0.201 0.182 0.185 -9.49% -8.15% 
FE^ 0.497 0.498 0.499 0.1996 0.4396 
RB 0.965 0.961 0.962 -0.37% -0.31% 
乳腺 0.987 0.99 0.991 0.32% 0.45% 
几 肉 .045 1.042 1.042 -0.24% -0.25% 
FAE .054 1.051 1.051 -0.30% -0.31% 
骨 松 质 200mg/cc .105 1.134 1.142 2.59% 3.38% 
骨 密 质 800mg/cc .406 1.425 1.433 1.37% 1.94% 
上 骨 密 质 1250mg/ce .62 1.607 1.602 -0.79% -1.13% 
3.2 与 SECT 单 参 数 比 较 
2.0 
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到 2. 拟 合 HU-SPR 曲线 


采用 120kVp 扫描 参数 ， 使 用 传统 SECT 单 参数 提取 方法 "， 拟 合 出 HU-SPR 曲线 (图 2)， 再 通过 


线 ， 计 算 组 织 等 效 模 体 的 SPR 值 与 误差 ， 并 于 表 5 中 双 参 数 提取 方法 计算 SPR 的 误差 作对 比 〈 表 6)。 
表 6. 与 传统 方法 误差 作 比 较 


DECT SECT 
组 织 等 效 材料 
DECT80/140 DECT100/140 SECT120 

吸 气 肺 -9.49% -8.15% -0.50% 
FEE 0.19% 0.43% 4.02% 

KA -0.37% -0.31% -2.49% 
乳腺 0.3296 0.4596 -0.4196 

几 肉 -0.24% -0.25% 1.72% 

FE -0.30% -0.31% -1.23% 
骨 松 质 200mg/cc 2.59% 3.38% 1.45% 
骨 密 质 800mg/cc 1.37% 1.94% 0.50% 


骨 密 质 1250mg/ce -0.79% -1.13% 0.19% 


mSECT120 


m DECT100/140 
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3. 两 种 方法 预测 的 SPR 偏差 对 比 
4 讨论 
质子 重 离子 放疗 作为 精准 放疗 手段 ， 近 年 来 , 引起 了 国内 外 放疗 界 的 广泛 关注 。CT 值 与 相对 阻止 本 领 
转换 曲线 对 于 质子 重 离子 放射 治疗 计划 制定 的 准确 程度 起 到 十 分 重要 的 作用 "。 本 研究 通过 采用 常规 定位 
CT 设置 不 同 管 电压 时 产生 的 不 同 能 量 的 X 射线 对 组 织 等 效 模 体 扫描 , 并 通过 高 低能 量 扫描 获得 的 两 组 CT 
值 ， 提 取 p。、Zerr 双 参 数 ， 并 通过 双 参 数 得 到 介质 的 相对 阻止 本 领 


本 文 方法 提取 的 参数 〈 表 4) pe 的 1 


偏差 .， 除 骨 松 质 (Bone 200mg/cc) 误差 达到 -8.39% 与 -9.70%2 


中 通 


r3 因为 多 孔 材 料 在 
偏差 的 | 


原因 ， 可 能 也 是 由 于 带 
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